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Abstract 

In this paper, an existence theory is established for ring-profiled optical vortex solitons. We 
consider such solitons in the context of an electromagnetic light wave propagating in a self-focusing 
nonlinear media and governed by a nonlinear Schrodinger type equation. A variational principle 
and constrained minimization approach is used to prove the existence of positive solutions for an 
undetermined wave propagation constant. We provide a series of explicit estimates related to the 
wave propagation constant, a prescribed energy flux, and vortex winding number. Further, on a 
Nehari manifold, the existence of positive solutions for a wide range of parameter values is proved. 
We also provide numerical analysis to illustrate the behavior of the soliton’s amplitude and wave 
propagation constant with respect to a prescribed energy flux and vortex winding number. 
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1. INTRODUCTION 


An exciting area of research in modern optics is the study of optical vortices. In a 


light wave, optical vortices are formed by wave dislocations or defects 


28l |. Its applications 


are found across numerous branches of nonlinear science, such as quantum information 


processing, wireless communications, and some not directly related to wave propaga 
instance condensed matter physics, particle interactions, and cosmology js-?. 
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38|, 1^. An interesting class of optical vortices are the ring-prohled optical vortices. Such 


vortices can be considered as a ring of light with a black spot at its center. In terms of a 
light beam, such black spots represent a zero light intensity. 

Of particular interest, is the theoretical description of a comp lex-valued light wave prop¬ 
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in a nonlinear media and governed by a nonlinear Schrodinger equation [l, 
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Rigorous mathematical treatments of such nonlinear problems 


present mat 

i bHid, I24, 


lematical challenges and have been considered by mathematical analysts J2 


33 


39 


40l |. Our interest is motivated by the work of Skryabin and Firth [36| 


and the mathematical analysis of Liu and Ren 


24j |. and Zhang and Yang 401. 


Consider the propagation of an electromagnetic wave in the longitudinal ^-direction over 
the transverse plane of coordinates {x, y) perpendicular to the z-axis. In dimensionless form, 
the evolution of the slowly varying electric held envelope, E, is modelled by the nonlinear 
Schrodinger equation |36l |. 


dE 


ylE + E{I)E = 0 , 


( 1 . 1 ) 


dz ' 2 

where is the Laplace operator over the transverse plane of coordinates. The function 
E depends on the total held intensity, I, i.e., I = and encapsulates the nonlinear 

properties of an optical medium. Examples for E{I) include 


1. E{I) = I (pure Kerr nonlinearity), 

2. E{I) = I — aP (cubic-quintic model), 

3. E{I) = /(I -|- al)~^ (saturable nonlinearity). 


where a G M is a parameter describing the nonlinearity saturation 12|]. The saturation 


constant may be dehned as a = Ini/ld, with Ini and Id denoting the nonlinear and dihraction 
lengths, respectively. Note that the saturable nonlinearity and cubic-quintic models reduce 
to the pure-Kerr nonlinearity when the saturation constant is zero. 
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We focus on spatial optical solitons 


0,0 


35 


37j . Spatially localized solutions of fll.ip . 


which do not change their intensity profile during propagation, can be described under the 
spatial soliton ansatz 

Eix, y, z) = u{x, y) exp(iK 2 ; + i(j){x, y)), (1.2) 


where u and 0 are real valued functions representing the soliton amplitude and phase, 
respectively, and k G M is the wave propagation constant. In view of fll.2p . the nonlinear 
Schrodinger equation fll.ip transforms into the coupled system 

I Vin - - 2[/s: - F{u^)]u = 0, 

[ Vi0 + 2V0-Vln(n) = 0. 

Under an appropriate balance between the nonlinear and diffraction lengths, the electro¬ 
magnetic radiation may become self-trapped and form a self-induced wave-guide. Solutions 
to fll.3p are referred to as self-trapped nondiffracting solutions in a self-focusing saturable 
nonlinearity. 

The existence theory in this paper seeks self-trapped positive radially symmetric solutions 
of fll.311 with a phase singularity at its center. Such solutions describe ring-profiled vortex 
solitons and can be found under the n-vortex ansatz 


u = u{r), (p = n9, r = 6 = arctan(?//a;). 


(1.4) 


where r, 9 are polar coordinates over and n G Z is the vortex winding number. 

Due to the presence of the vortex core or, equivalently, the regularity of u at r = 0, we 
impose the condition m( 0) = 0. Moreover, such ring-like beams remain localized. Thus allow¬ 
ing us to mathematically impose the “boundary” condition u{R) = 0 for i? > 0 sufficiently 
large, where R represents the distance from the vortex core. 

Using (II.4p . and in a saturable nonlinear media, the system (II.3p reduces to the n-vortex 
equation 


2ru^ 


tR 

{rUr)r - U + — -r 

r 1 -|- au^ 


2Kru = 0 , 


m(0) = 0, u{R) = 0. 


(1.5) 


An important parameter characterization of spatial solitons is its energy flux. Using the 
n-vortex ansatz, the soliton energy flux is defined as 

pR nR 

Q{u) = / \E\‘^rdrd9 = 2% ru^dr. (1.6) 

Jo Jo 
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The rest of the paper is summarized as follows. In section 2, we give a necessary condition 
for the existence of nontrivial solutions. In section 3, we treat fll.Sp as a nonlinear eigenvalue 
problem and use a constrained minimization approach, subject to a prescribed energy flux 
constraint, to prove the existence of positive solution pairs {u, k). In section 4, we prove the 


existence of positive solutions for a wide range of parameter values over a Nehari manifold. 
In section 5, we supplement our results by using a finite element formalism to compute the 


soliton’s amplitude and wave propagation constant for a prescribed energy flux. A summary 
is provided in section 6. 

2. NECESSARY CONDITION FOR NONTRIVIAL SOLUTIONS 

The n-vortex equation (11.hh may be viewed as the Euler-Lagrange equations of a cor¬ 
responding action functional. For a sufficiently large distance R, we prove that the action 
functional is indehnite, and, therefore, a direct minimization approach is not possible. 
Consider the action functional —)■ M defined as 


with \n\ > 1 and a > 0. hf is the completion of 



X={ue C^[0,i?]|M(0) = 0 = u{R)} 


( 2 . 2 ) 


(the space of differentiable functions over [0, R] which vanish at the two endpoints of the 
interval) and is equipped with the inner product 



(2.3) 


We may treat H as an embedded subspace of composed of radially symmetric 

functions enjoying the property m( 0) = 0 for any u E H, where Br := {(a:, y) E < R^}. 


We are interested in finite energy solutions of (11.511 . From the inequality 


ln{l + x^) < x'^ for all a: G M 


(2.4) 


and 



(2.5) 
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we note that the norm induced by the inner product fl2.3p on H guarantees that all terms 
in the functional stay hnite. In other words, there is a constant C > 0, such that 
< C'||m|||/- For convenience, we dehne the ‘energy’ functional as 

1 rR 


£{u) = 


rul + -u^ + r ln(l + au^) \ dr. 


( 2 . 6 ) 


Theorem 2.1 Ifu is a nontrivial finite energy {8{u) < cxo) solution of the n-vortex equation 
fll.Sp . then the wave propagation constant must satisfy 


_i rl + 
K < a - 


2 i ?2 ’ 


(2.7) 


where ro(~ 2.404825) is the first positive zero of the Bessel function Jo JSfiJ . 


Proof. Suppose if liminh^o {ru{r)\ur{r)\} 7 ^ 0, then there is an e > 0 and ro G (0, R] such 
that ru{r)\ur{r)\ > e for all r G (0,ro). However, 

,2 \ 1/2 / I'TQ \ 1/2 


00 = 


'0 


ro , rro 

-dr< / 

r Jo 


u\u 


r\dr < ij — dr\ ij ru^drj , ( 2 . 8 ) 


which contradicts the hnite energy condition. Hence, liminf {ru{r)\ur{r)\} = 0 and we can 


r—^O 


extract a sequence {r^} such that —)■ 0 as j —)■ cx) and 

liminf {rju{rj)\ur{rj)\} = 0 . 

r —>-0 


(2.9) 


Multiplying fll.Sp by u, we integrate by parts over the interval [rj,R] and let j —)■ cx). 
Appealing to fl2.9p . we have 

2 , i 2 2 2 ru^ . 

ru dr = < —u + 2nru - 7 / dr. 

Jo Jo { r 1 + au^ 

Using F{u) = i < a~^ and fl2.5p . we get 


( 2 . 10 ) 




ruf.dr > —2 { a 


-1 


n 


rR 


/ ru^dr. 
2i?2 ) ./o 


( 2 . 11 ) 


Treating u as a radially symmetric function in and using the Poincare inequality 

rR 


r2 rR 

ru^dr < / ruidr, 

n Jo 


( 2 . 12 ) 


5 
















with ro as defined in fl2.7l) . gives 


0>-2('a-'-!^-d TrAr. 




Let 


_i 'n? + 

a ■.= a -—- K. 


2 i ?2 


(2.13) 


(2.14) 


If (j < 0, then (I2.13P gives a contradiction. Therefore, we must have a > 0. □ 

As a consequence of the following lemma, when the distance from the vortex core R is 
sufficiently large, we prove that the functional fl2.ip is indefinite. 


Lemma 2.2 Let k < a ^ and fd > 0. If the distance from the vortex core satisfies 


R > 


12(1 + n2(21n2 - 1)) 


a ^ — K 


then there exists an element Uq & H such that 

\\uo\\h >/3 and X«,(mo) < 0 . 


(2.15) 


(2.16) 


Proof. Set R = 2a and define 

uo{r) = 

By direct calculation we obtain. 


b 

a ’ 


0 < r < a. 


-(2a — r), a < r < 2a. 


(2.17) 


f*2a 


ruidr = 


/o 

f*2a 


rul^^dr = 26^, 


'*2a 


-u^dr = 26^ (2 In 2 — 1), 


(2.18a) 

(2.18b) 

(2.18c) 


[ r \n{l + aul)dr = 2a'^ (lrL{l + ab'^) — 2 -\—^tan ^(\/a&)l , (2.18d) 

Jo V \/«o / 

where Mo,r := (wo)r- Similar to Lemma 3.3 in |39[, we note that uq is obtained as the limit 
of a Cauchy sequence in H, consequently, Uq belongs in H. Using (2.19), we get 


||no||?, = 46^1n(2), 

Xk(mo) = b ‘^ “ 1) ~ ~ 


(2.19) 

( 2 . 20 ) 


2a^ 


ln(l + q;6 ^) — 2 H-p=-tan ^i^/ab) 

^Jab 


6 














with a = R/2. For any e > 0, choose b sufficiently large such that 




ln(l + — 2 H—^tan ^(^/ab) 

^^ab 


< e. 


Hence, 


'^k{uq) < 6^ ( 1 + n^(2 In2 - 1) + ( e - ^{a ^ - k) ] 


( 2 . 21 ) 


( 2 . 22 ) 


In order for the right hand term of (12.221) to be negative, let e = -(a ^ — k) and R be such 

3 

that 


R > 


/12(1+ n2(21n2 - l))y 




a ^ — K 


(2.23) 


For any j3 > Q, choose b satisfying fl2.2ip and, such that, ||mo||// > /d- With these values of 
b and R, we have Xk{uo) <0. □ 

From inequality (I2.22jl . if the distance from the vortex core is sufficiently large, then 
^k(wo) —oo as 6 —)• cx) for K < Therefore, for a sufficiently large distance R and 

K < a~^, the functional is indehnite and as such, a direct minimization is not possible. 

3. EXISTENCE VIA CONSTRAINED MINIMIZATION 


Using a variational principle and constrained minimization problem, we prove the exis¬ 
tence of positive solutions of the n-vortex equation (II.hh . In this scenario, the wave prop¬ 
agation constant k is undetermined and appears as a Lagrange multiplier. We provide a 
series of explicit estimates for the wave propagation constant, vortex winding number, and 
a prescribed energy flux. 

We view (11.5p as a nonlinear eigenvalue problem 

2ru^ 


n 


{rur)r - u -h . „ 

r 1 -|- au^ 


= 2Kru, 


m(0) = 0, u{R) = 0. 


Dehne the action functional X as. 


I ( r)2 


ru^ H- w — 2a + 2a ln(l -f au^) )■ dr, 


and the soliton energy flux constraint functional Q as 

.R 

Q{u) = 271 ru^dr. 


(3.1) 


(3.2) 


(3.3) 
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Consider the nonempty admissible class 


A = {u(r) is absolutely continuous over [0,i?],M(0) = u(R) = 0,£^(u) < cxd} , (3.4) 

where S^(u) is a dehned by (I2.6p . In order to prove the existence of a solution pair (u, k), it 
suffices to show that a solution to the following exist: 

7] = inf {I(u) : Q(u) = Qo > 0}, (3.5) 

u&A 

where Qo is a prescribed value for the energy flux and k is the Lagrange multiplier. 

Theorem 3.1 Consider the n-vortex equation fll.5p . describing ring-profile vortex solitons 
in a self-focusing saturable nonlinear media, subject to the prescribed energy flux Q{u) = 
Qo > 0 and finite-energy condition £{u) < oo, defined by fl3.3p and fl2.6p . respectively, with 
parameters \n\ > 1, a > 0, and R> 0. 

(i) There exists a solution pair {u, n) satisfying u{r) > 0 for r G (0, R). 

(a) //n^ + 2r^K > 0 forr G [0, R], then there exists no nontrivial small-energy-flux solution 
satisfying Q{u) = Qo < 1/4. 

Proof, (i) From the prescribed energy flux, it follows that the functional R{u) satishes 

R R ^ ^ 

>^[ ruldr + / —dr - ^^Qo- (3.6) 

2 J Jr 27r 

0 0 

As a result, the minimization problem (13.5p is well-dehned. Let {uj} be a minimizing 
sequence of fl3.5p . i.e., choose a sequence of functions {uj} in A such that 

T{uj) —^7] as j ^ oo and X(mi) > X{u 2 ) > ■ ■ ■ > rj. 

Since {uj} minimizes (13.5p and using (13.6p . there exist C > 0 independent of j such that 

rR i-R y2 

C> ru^- dr / —dr, (3.7) 

Jo Jo r 

where Uj^r ■= 

The distributional derivative of u satishes Hwlrl < \ur\, and the functionals X and Q are 
both even, i.e., X(m) = X(|w|) and Q{u) = Q{\u\). As a consequence, we assume that the 
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sequence {uj} consists of nonnegative valued functions. Moreover, we take these functions 
to be radially symmetric over the disk Br and vanishing on its boundary. 

From fl3.7p and fl2.5p it can be seen that the functions Uj belong in under the 

radially symmetric reduced norm, 


pR rR 

||m|P := / rv?dr + / ruldr. 

Jo Jo 

Using 03.71) and 02.Sp . we show the sequence {uj} is bounded in WQ’‘^{Bii), 

pR pR pR y2 pR 

ll'*^ 7 ll^= / ru^.dr + / rvJ„dr<R^ / -J-dr + / + 1)U. 

Jo ^ Jo Jo r Jo 

Without loss of generality, since we are in a reflexive space, we may assume the weak 
convergence of {uj} to an element u G Wo’‘^{Br). As a result, it now suffices to show that 
Uj converges to a minimizer of 03.5|) and belongs in A. 

From the compact embedding W^’‘^{Bfi) CC Lp^Br) for p > 1, Uj ^ u strongly in 
U’{Bpi) as j —?• oo. Hence, u is radially symmetric and satisfies the boundary condition 
u{R) = 0. 


In view of 03.71) and using Fatou’s lemma, we get 


r-H 


r-H 


ru^dr < liminf / ru^^dr, 

Jo 

U^j 

/ —dr < liminf / —dr, 


'0 


r 


J^OO 


>0 


r 


i-R 


i-R 


rln(l + cra^)dr < liminf / rln(l + aMHdr, 

j^oo Jo 


(3.8a) 

(3.8b) 

(3.8c) 


where the finiteness of the right hand side of 03.8cp follows from 02.4p and 02.Sp . Therefore, 
from 02.sp and (3.8), we get the weak lower semi-continuity of the functional X, i.e.. 


X(n) < liminfX(nj 




(3.9) 


Using 03.9p . together with 03.Sp . gives X(m) = p. Further, note that Q{u) = lim Q{uj) = Qo- 


J^OO 


Moreover, the finite-energy condition also holds from 03.7p and 03.8ap - 03.8cp . In particular, 
ru^, rul, and ln(l -|- au^) are all in X(0, R). 

To show that m( 0) = 0, we follow as in [39|. Let {uj} be a sequence in IU^’^(e, R) where 
e G (0,72). For any e G (0,72), {uj} is bounded in IU^’^(e,72). The compact embedding 
IU^’^(e, 72) CC C[e,72], then gives Uj —)■ u uniformly over [e, 72] as j oo. Thus, for any 
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pair ri, r 2 G (0, R) such that ri < r 2 and using C from fl3.7l) . we get 


\u%r2)-u%n) \ = 


n2 


{u]{r))rdr 


fri 

pr2 


< / 2\uj{r)uj^r{.'r')\dr 

J ri 

/ rr2 \ 1/2 

< 2 ( / ru^^^{r)dr j 

V^ri / \^ri 


-ra \ 1/2 


M^(r) 


dr 


(3.10) 


Since Uj —)■ n uniformly as j —)■ cx3, we take j —)■ cxo above to get 





(3.11) 


The right hand side of (13.lip goes to zero as ri,r 2 —)■ 0; since ^ is in L(0,i?). Hence, the 
following limit exists, 

^0 = limn^(r) = 0. (3.12) 

r-^-O 

As a consequence, the boundary condition m( 0) = 0 is achieved. 

Therefore, the function u, obtained as the limit of the minimizing sequence {uj}, is a 
solution to the constrained minimization problem fl3.5p . and there is a real number k such 
that {u, k) satishes (13.Ih . 

Further, we may suppose that there is a point tq G (0, R) such that u{ro) = 0. Since 
tq would be a minimum point for u{r), we have Mr-(i"o) = 0. However, by the uniqueness 
theorem of the initial value problem of ordinary differential equations, u{r) = 0 for all 
r G (0,i?), thus contradicting the energy flux constraint Q{u) = Qq > 0. Hence, u{r) > 0 
for all r G (0,i?). A standard bootstrap method may then be used to conclude that u is a 
classical solution of (II.5p . 

(//) Let {u, k) be the solution pair obtained in part (/). Using (I2.10p . we 

get 




ruzdr > 


r-R /^2 


-R 


— + 2 r I ru^dr 


ru^dr. 


(3.13) 


We treat m as a radially symmetric function dehned over with its support contained in 
the disk Br. From the classical Gagliardo-Nirenberg inequality over R^, we write 


r-R 


fR 


-R 


ru^dr < dvr / ru^dr / rutdr. 


(3.14) 
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As a result, we get 


.R 

— / ruldr > 


r-H /^2 


fR 


-R 


+ 2k ru^dr — Sir ru^dr / rutdr. 


(3.15) 


Rearranging the terms in fl3.15p . and using the prescribed energy flux constraint gives 

rR 


(4(5o - 1) / ruldr 


/^2 


+ 2k ru^dr > 0. 


Hence, u = 0, if 


1 72^ 

Qq < - and — + 2 r > 0 for r G (0, R], 


(3.16) 


(3.17) 


as claimed. 


□ 


Theorem 3.2 Let (u, k) he the solution pair of the n-vortex equation fll.Sp obtained in 
Theorem 3.1 with k as the wave propagation constant. 


(i) The wave propagation constant satisfies 


K >a 


-1 


Jj(l + n=(21n2-l)) 


tiR^ 


"" V ^ J + 


SaQo 


tan 


-1 


(a) If the vortex winding number satisfies |n| > Qo/n, then k <0. 


SaQo 

7rR2 


(3.18) 


(Hi) For K > 0, the solution pair (u, k) satisfies 


u 


< CKexp(-V^r), 


for r sufficiently large and > 0 is a constant depending on k only. 


Proof, (i) To obtain a lower bound for k, we rearrange fl2.10p and write 


rR 

K / ru^dr = — - 


-R 


2 2\ 1 . 
ru„ H- u ar + / - -dr. 


'0 


1 + au"^ 


(3.19) 


Choose uq G a, satisfying Q(uo) = Qo > 0. Since w is a solution to the constrained 
minimization problem fl3.5p . whose existence was proved in Theorem 3.1, we have T(u) < 
T(uo). As such, we get 


-R 


n 


-R 

ru) H-) dr < X{uo) + a~^ / ru^dr — a~'^ I rln(l + au^)dr. 


-R 
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Inserting the above into fl3.19p and using Q{u) = Qo > 0, gives 


_ 1 Stt - > _2 27r 

K>-a - —X(mo) + a — 
Vo Vo 


r-H 


r ln(l + av?)dr + 


27r 


r-i? 


rti 


Qo Jo lx 


-dr. (3.20) 


Using the inequality 


ln(l + x) > 


X 


1 + a; 


for X > 0 , 


(3.21) 


fl3.2ip may be rewritten as 


27r 


2ti V? + au 


K > —a — —X(mo) + a — / r 


27r 


Qo 


■dr = ——X(mo). 


Qo Jo 1 + ”” Qo" 


(3.22) 


From fl2.17p and (2.19), we obtain 

X(mo) = &^ + In 2 — 1 )- a~^a^b 


2 


(3.23) 


+ 


2 a^ 




ln(l + ab'^) — 2 + 


2 tan ^{^/ab) 

^/ab 


Using (I2.18al) and Q(mo) = Qo, we get 6 ^ = 3Qo/(vri?^). Inserting (I3.23p in (I3.22p . we arrive 
at 


6 


K 


>a-^ - — {1 + n\2\n2 - 1)) 


(3.24) 




In (l + ab‘^) — 2 H—p=- tan ^ (^/ab) 
^ yjab ^ 


which is the desired lower bound. 

[ii] Let (m, k) be the solution pair obtained in Theorem 3.1. Using Schwartz’s inequality, 
and m( 0 ) = 0 , we get 


u 


Jr) = / 2u{p)up{p)dp < 2 / pu {p)dp 


r-R 


1/2 / pR 


U^ip) 


lo P 


dp 


1/2 


(3.25) 


Multiplying fl3.25p by ru^, integrating from 0 to R, and using the constraint Q(m) = Qo > 0, 
we obtain 


-i? 


ru'^dr <^([ pu^p{p)dp 


TT 


1/2 / pR 


u\p) 


lo P 


dp 


1/2 


(3.26) 


Using inequality, ab < ea^ + X for every a, 6 G M and e > 0, gives 


r-R 


-R 


ru^dr <e pui{p)dp + 


Q^ f^u\p) 


e47r2 Jo p 


dp. 


(3.27) 
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From fl3.27p and in fl3.19p . we get 

rR 


l+au^ 

-R 


K 


ru^dr < 


'0 




- , Ql 


u\ 
— dr. 

r 


(3.28) 


'0 


We choose e = | in fl3.28p and conclude that k < 0 whenever |?7,| > QQ/n. 

{Hi) The exponential decay estimate follows from an application of the maximum principle 
and a suitable exponential comparison function. We rewrite (13.11) as 


Am = Urr H- Ur = 


n 


2u^ 


1 + au'^ 


+ 2k\ u. 


(3.29) 


It then follows, 

> 2mAm = 2 f ^- + 2i^u^ >a(k - - —wl (3.30) 

\r^ 1 + aw"' / \ 1 + au^J 

By the continuity of w on [0, R] and the boundary condition u{R) = 0, for any e > 0 there 
is an > 0 such that 


Am^ > 4 (k — e) for every r G [R^, R\. 
Define the comparison function ^ : [0, i?] —?• M as 

^{r)=Ce-^^, C,a>0. 


(3.31) 


(3.32) 


Hence, 

A^ = - —. (3.33) 

r 

Subtracting fl3.3ip and fl3.33p . for every r G [R^, R], 

A(m^ — ^) > 4 (k — e) > 4 (k — e) u^ — (3.34) 

For any k > 0, we choose e satisfying 0 < e < k and = 4 (k — e), to get 

A(m^ — 0 ^ a^(M^ — C) for every r G [i?e, R\. (3.35) 

Let C in fl3.32p large enough so that ^ < 0 for r = Denote C by to emphasize 

its dependence on e. Since, ^ 0 as r —)■ R~, and applying the maximum principle, we 

conclude that < 0 for all r G [i?^, R]. For simplicity, let e = m/ 2 to obtain 

< CKexp(—\/^r) for every r G [Rk,R], (3.36) 


13 




















where > 0, and depends only on k > 0. □ 


Remar 

at infinity 


cs. Beam confinement reqnires the exponential decay of the soliton amplitnde u 


36| . From the exponential decay estimate given in Theorem 3.2 {Hi), we see that 


this occnrs for k > 0. Theorem 3.2 {ii) states that the vortex winding nnmber mnst satisfy 
\n\ < QqItt. On the other hand, from Theorem 3.1 {ii), if the wave propagation constant is 
positive, then the prescribed energy finx mnst satisfy Qo > As will be seen in Section 6, 
the condition on the prescribed energy finx, Qo > is not snfficient to conclnde that the 
propagation constant is positive. 

When the prescribed energy finx Qq is fixed, and the distance from the vortex core R goes 
to infinity, the necessary condition given in Theorem 2.1, together with the lower bonnd for 
the wave propagation constant, gives the ineqnality 


0 < K < a 

which is in agreement with the resnlts of Skryabin and Firth 
solntions of their model. 


(3.37) 


36l |. for any self-trapped 


4. SOLUTIONS ON THE NEHARI MANIFOLD 

Recall the action fnnctional : FT —>■ M defined as, 

|rM^ -I- —u^ — 2{a~^ — K,)ru^ + 2a~‘^r ln(l -|- au^) | dr, (4.1) 

where \n\ > 1 and a > 0. 

Standard argnments show that X^ G C^(iF, M). Also, 

f TUU 1 

{I'{u),u) = < rUrUr H- uu — 2{a~^ — n)ruu + 2a~^ -^ > dr, Wu G H, (4.2) 

Jo [ r 1 au^ J 

and (•,•) denotes the nsnal dnality between ff and its dnal space iF“h Let : iF —)■ M be 
defined by 

7 k{u) = l(I'Ju),u} = ^ f I ru^ + —u^ - 2{a~^ - K)ru^ + 2a~^ | dr. (4.3) 

2 2 Jq I r 1 -h au^ J 

For a fixed propagation constant k, define the Nehari manifold A4 as 

M = {ue H\{0} : 7,(n) = 0}. (4.4) 
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If M G i/\{0} is a critical point of X^, then for every v G H, {Xk{u),v) = 0. Setting v = u, 
it follows that 7 k(m) = |(Xk(m),m) = 0. Hence, u G A4, and the Nehari manifold contains 
all nontrivial critical points of X^ on H. 


Lemma 4.1 u E H is a nontrivial critical point ofl^. if ond only if u E M. and is a critical 
point ofX^lM- 


Proof. The forward implication follows directly from the dehnition of the Nehari manifold. 
We now justify the other direction. For every u E A4, 


^ 2 2 -1 \ 2} j o -1 f 

ruz H- — 2{a ^ — K)rvf > dr = —2a ^ / <- - 

r j Jo U + «« 


dr. 


(4.5) 


By definition of 7 ^, for every u E A4, we get 


(7k(w),«) = 




./o 
= 2a 


ey Th ey . _ 1 , 9 _ 1 VU 

ruz. H- u - 2(a - uVu + 2a — -^ 

r (1 + au^Y 

r{^( 1 -l|)dr<0. 

0 t 1 + ecu VI + CtM 


dr 


(4.6) 


For any critical point Uo G Ai of X^lx, there exists a Lagrange multiplier, G M, such that 
{X'Yuo),u) = ^( 7 ^(mo), F), for every u E H. Using the definition of Al, gives 


0 = 7»(mo) = = ((7Ut‘a),'Ua}. 


(4.7) 


As a result, using fl4.6|) . it follows that ^ = 0. Therefore, the critical points of XY\m ^^ee also 
the critical points of X^. □ 

Lemma 4.1 indicates that A1 is a natural constraint for X^. From the necessary condition 
obtained in Theorem 2.1, the Nehari manifold contains no critical points of X^, when k > 
a~ - 2 ^^- Hence, we consider the case n < a~ - 


Lemma 4.2 If the distance from the vortex core satisfies 

'6(1+ u2(21n(2) - 1)) 


a 


-1 


K 


< R, 


(4.8) 


then the Nehari manifold is not empty. 


Proof. Define 

F(f,M) 


r-R 


n 


ruf. H- — 2{a — K)rm + 2a 


, -1 


ru 


1 + aRu"^ 


dr. 


(4.9) 


15 












As a result, 'ynitu) = t‘^T{t,u). For any u 7 ^ 0 and k > —, we get 


rR 


n 


r(0, u) = - / < rul H -+ 2Kru^ > dr > 


rg + 
2i?2 


+ K 


fi? 


ru^dr > 0. 


(4.10) 


Note that 


-R 


t^QO 


r(oo, u) = lim r(t, n) = - / <( ru^ + —u^ — 2{a ^ — K)ru^ \ dr. 


(4.11) 


Substituting uq as defined in Lemma 2.2, in fl4.1ip . we get 

1 


r(cxo, uq) = h [1 + n (2 ln(2) — 1)- (a — k)R 

6 


(4.12) 


Selecting R as in fl4.8p . we get r(cxo,Mo) < 0. Hence, there exists a to > 0 such that 
r(to,Mo) = 0 and, it follows that, 7K(to'Wo) = = 0. Therefore, to'^^o is in fhe Nehari 

manifold A4. □ 


r 7*^ I 

Note that, for every k > — 




(7"(0 ), m) = / <( rul + —u^ + 2Kru‘^ '> dr > 2 


rg + 
2i?2 




+ K / ru^dr > 0. (4.13) 


Hence, u = 0 is a strict local minimum of 7^ and, as a result, an isolated point in Ad U {0}. 
Thus, 0 ^ dM.. Therefore, for all u G Ad, there exists a constant Ci > 0, independent of u, 
such that 


\u\\h > Cl- 


(4.14) 


Lemma 4.3 There exists a constant Co > 0, such that 


( 7 (.(m),m) < —C 2 < 0 for all u G Ad. 


(4.15) 


Proof. Let be a sequence in Ad such that 


fR 


lim / ^ „ 

j^Qo /g 1 + au 


^—^dr = 0 . 


(4.16) 


Hence, {uj}fLi is either bounded or unbounded in H. If is unbounded, then there 

is a subsequence such that pk := —)• 00 as A; —)■ 00 . Let Vk = and hence 
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= 1 - If necessary, passing to a subsequence, Vk v m H, Vk ^ v in and 

nfc(r) —)■ v{r) a.e. r G [0, i?] as /c —)■ oo. Consequently, 

r 

k^oo 


TV?‘- 

0 = lim / ^= lim 


rvl 


lo 2 


dr = -a R > 0 , 


(4.17) 


fc^oo Jq 1 + aUj^ 

which is a contradiction. Hence, must be a bounded sequence in H. If necessary, 

passing to a subsequence, we have Uj u in H, Uj —)■ u in and Uj{r) — )■ u{r) a.e. 

r G [0, i?] as j —)■ cxo. Using the dominated convergence theorem, we get 


0 = lim / 

j^oo Jq 1 + au 


^ ru] ru^ 

dr = ' 


Iq 1 “t“ (yu' 


;dr, 


(4.18) 


which gives u{r) = 0 a.e. r G [0, R]. Since is in A^, we get 


»R 


lim / 

j^oo Jq 


\ dr = lim <| (a ^ — k) / ru^-dr — 2 a ^ ‘ ^ 


j-)-oo 


'0 


'0 


1 + av?^ 


dr 


= 0 . 


(4.19) 


It follows for any \n\ > 1, 

0 = lim [ dr > lim > 0. (4.20) 

J^ooJq r i^oo" 

Hence, Uj —)■ 0 in id, which contradicts fl4.14p . Therefore, there exists a constant C 3 > 0 
such that 



2 

ru 

1 + au‘^ 


dr > C 3 


for every u G A4. 


Using (14.op . and Holder’s inequality, we get 


(7k('^),u) 



ru 


(1 + au 


2\2 


dr 




2 

ru 

1 + au'^ 



2 



—C*2 < 0. 


(4.21) 


(4.22) 


□ 


Lemma 4.4 The set A4 is a paracompact and complete topological space. 
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roof. The paracompactness of the set M follows identically to the proof in Lemma 3.6, 


24| . To show that A4 is complete, we let be a sequence in A4 such that uj —)■ u 


m 


H as j ^ oo. From the compact embedding, CC Lp^Br) for p > 1, we note that 

Uj —)■ u strongly in Lp{Br) as j —)■ oo. Hence, 


r-iJ 


j->oo 


0 = lim 'yniuj) = ^ I “i — 2(a ^ — K)ru^ + 2a 




ru 


1 + au"^ 


dr 


(4.23) 


Since I \u\\h > C 


= 1 k{u). 

> 0, by fl4.14|) . we conclude that u G A4. □ 


Similarly to 2^, using Lemma 4.3, we may deduce that A4 is a regular C^-Banach 


manifold and, moreover, using Lemma 4.4 that is a Finsler manifold. We now look for 
nontrivial solutions of the n-vortex equation 01.51) . as critical points of restricted to the 
manifold A4. 


Lemma 4.5 Let {uj}JLi be a sequence in M. such that {Xk(mj)} is bounded. Then the 
sequence is bounded in H. 


Proof. Let be a sequence in M. such that {Xk(mj)} is bounded. Hence, 

0 = - i^)ru] + 2a"^:—, 

2 Jo I r 1 + otuj) 

and there exists a constant /? > 0 independent of j such that 


(4.24) 


\'^K ) I 


fH 


Tl 

H-— 2(a“^ — k)™^ + 2 a~'^r ln(l + aw^) [> dr 


Using 04.24P and 04.25p . we get 

rR 


ln(l + auy ^ 


^ 1 + au] 


rdr < 0 ^( 5 . 


< f3. (4.25) 


(4.26) 


JO ^ 

Assume the sequence {uj}^^ is unbounded in H. Let pj = Then, pj oo as 

j —)■ oo. Let Vj = —. Then ||r',||ir = 1. Hence, without loss of generality, we suppose that 

J Pj M J M 


Vj V 


(021, 


in H, Vj —)■ n in Lp{Br) for every p > 1, and Vj{r) v{r) a.e. r G [0,i?]. From 


1 — II'^iIIh ^ 


< 


^ 2 2 I 7 

rVj^^ + —Vj } dr = 


.-1 


2(a — K)rVj — 2 a „ 

01 ^ l + aw? 


dr 


-i? 


{2(a ^-K)r'i;|}dr = 2(a ^ - R)\\vj\\l^B^y 


(4.27) 
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Letting j —>■ cx) in fl4.27D . we get 1 < 2{a ^ Hence, v ^0 a.e. r G [0,R]. Let 

H = {r e [0,R] : v{r) ^ 0}. Then |H| ^ 0. Using fl3.21l) . it follows 


r-H 


ln(l + au^) — - -^ 

^ l + QfM? 


rdr > < ln(l + au^A 


au] 


1 + auj 


rdr. 


(4.28) 


Note that as a result Vj{r) = ^ ^ g a.e. r e as j —)■ cx3. Hence, \uj{r)\ — )■ cxo 


CXU '^ 

and ln(l + aUj(r)A - - ——>■ oo a.e. r G H. Applying Fatou’s lemma and fl4.26p . 

1 + 

we get the contradiction 




a‘^(3 > liminf / < ln(l + auA 


J^OO 


au. 


1 + av?j 


> / liminf <( ln(l + atiU — ^ 


1 + avid- 


In 

Therefore, the sequence {uj}^^ is bounded in H. □ 


rdr 

rdr = oo. 


(4.29) 


Lemma 4.6 satisfies the Palais-Smale condition onAi, namely, if is a sequence 

in A4 such that {X^iuj)} is bounded andXfi\'j^{uj) —)■ 0, then there exists au E M. such that 
Uj —)• u (strongly) in H. Moreover, u is a critical point ofXfi\M- 

Proof. Suppose that {X,.,_{uj)} is bounded. Then, from Lemma 4.5, is bounded in 

H. Without loss of generality, there is a sequence {uj}°Ai such that Uj u in H, Uj u 
in L^{Br) for every p > 1, and Uj{r) — )■ u{r) a.e. r G [0, i?]. For every v E H, we have 


{Kiuj),v) = 


i-R 


r-i? 


^ _1 ^ rUjV 

rUirVr-\ - UnV — 2{a — njrUiV + 2a - - —^ 

r ^ ^ ^ ^ l + au] 


dr 




n / \ _i ruv 

rUrVr H- uv — 2(a — njruv + 2a - - 

r 1 + 


dr 




(4.30) 


and 


-i? 


n 


(7K('^i) w) = / "I ruj^rV H- UjV — 2(a — n)rujV + 2a 


rUjV 






n 


.-1 


(1 + auA"^ 


ruv 


rUrVr H- uv — 2(a — K)ruv + 2a 

r (1 + ati"*)^ 


dr 

dr 


= (7l(«)w)- 

Using the dehnition of X^K^f, there exists a sequence in R such that 

^0 in as j -E oo. 


(4.31) 


(4.32) 
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If { 0)^1 is unbounded, then there exists a renamed subsequence such that ^ oo. 

From (| 4.32 p, 0(7 kK),w) ^ and, hence, {'y'^{uj),u) 0. As a result, {'y'^{u),u) = 

0. However, using fld.lSp . we get 


0 > -C 2 > lim (7«(Mi),Mi) > = 0 , 


J-fOO 


(4.33) 


i=i’ 


a contradiction. Therefore, is bounded and contains a renamed subsequence 

such that ^ ^ as j —)■ cxo. Consequently, fl4.32p implies 


Kiu) - = 0 - 


Suppose ^ = 1. From fl4.34l) . 


0 = {X'^u) - -f'^{u),u) = 2a 




-1 



ru 


Iq ( 1 + 

which implies that u{r) = 0 a.e. r G [0,i?]. Thus, 


0 = IniUj) = -{X'^{Uj),Uj) 


1 - 


1 + au‘^ 


dr, 


(4.34) 


(4.35) 


(4.36) 





9 n 0 , ru, 

rUj ,,. H- Uj — 2{a — K)ruj +2a ^ 


1 + au'j 


dr, 


'0 t r 

and using the dominated convergence theorem, we conclude that 0 , which is a 

contradiction to fl4.14p . Therefore, 7 ^ 1. 

The boundedness of {uj}'^-^ and fl4.32p . gives 


f-i I \ uj — "^1 |iy —t 0. 


(4.37) 


\{K(Uj) - -U)\< \ \X'^{Uj) - 

then applying (I4.34p . we have 

(Kiuj) - - u) =(X1(mj) - - x'^iu) + ii^{u),Uj - u) 

= {KiUj) - X'^{u),Uj -u)- - ii^{u),Uj - u) 

= {'I-'KiUj) - X'^{u),Uj -u)- - i^{u),Uj - u) 

- (Ci - 0(7 kH, - u). (4.38) 

Using the dehnition of and 7 ^, we obtain 


R 


n 


{X'^{uj) — X'^{u),Uj — u) = J ^r{uj^r — Ur) H- {uj — uY — 2{a — K)r{uj — u)^ dr 


+ 2Q! 


-1 



Ui 


u 


'0 


1 + au^j 1 + q ; m | 


{uj — m) rdr 


(4.39) 
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and 


r ^2 


i'y^iuj) - 7 ^( m ), - u) = J - Mr)^ + —{uj “ “ 2(q! ^ - K)r{uj - uY \ dr 

rR 


+ 2 a 


-1 


Ui 


u 


Iq 1 V (1 + au^jY (1 + au^jY 

Applying the dominated convergence theorem, we get 


{uj — m) [■ rdr. (4.40) 


f-R 


Ui 


u 


0 t V1 + 


(uj — m)|- rdr —)■ 0 


(4.41) 


and 


fi? 


Mo 


u 


{uj — m) > rdr 0. 


(4.42) 


Jo + Y + c^u'jY, 

Since Uj ^ u m H, Uj —)■ w in L^{BY) for every p > 1, and Uj{r) u{r) a.e. r G [0,i?], 
from eqnations 04.371) - 04.421) and 02.Sp . we conclnde 


i-R 


(1 —.^) / r(Mj^r — Wr)^dr —)■ 0. 


(4.43) 


Therefore, \\uj — M||iy —?• 0 in H. From the completeness of A4, Lemma 4.4, u G Ai. Then, 
04.34P implies that w is a critical point of T^m- D 


Theorem 4.7 Let the distance from the vortex core satisfy 04.81) . For each propagation 
constant in the interval 


+ Vq _o -I- Tq 
a 


(4.44) 


2 i ?2 ’ 2/22 

there exists a solution pair (w, k), satisfying u{r) > 0 for r G (0, R), to the n-vortex eguation 

jia. 


Proof. Let u G M.. Hence 7 k (w) = 0 and 




n 


i-R 


ruj H- > dr = < 2(q! ^ — K)ru^ — 2a 


.-1 


ru 


1 -|- q;m2 


dr. 


(4.45) 


Inserting 04.45P into Xk and using 03.2ip . gives 


-i? 


iJu) = a 


-2 


an 


ln(l + au)- 

'0 i 1 + au^ 


rdr >0 on A4. 


(4.46) 


Thus the functional Xk is bounded below on Ad. As a consequence of Ekeland’s variatonal 


principle 


IJ, ll8j and Lemma 4.6, there exists a w G A4 such that Xk(m) = inf{XK(M)|M G Ad} 
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and Xk|^(m) = 0. By Lemma 4.1, u is also a critical point of and, therefore, a solntion 
of the n-vortex eqnation fll.51) . 

We nse the evenness of the fnnctional X^ to get a positive solntion. Moreover, u{r) > 0 
for all r G (0,i?). Snppose there is a point tq G (0, i?) snch that u{rQ) = 0. Then tq wonld 
be a minimnm point for u{r) and = 0- By the nniqneness theorem of the initial valne 

problem of ordinary differential eqnations, u{r) = 0 for all r G (0, R), thns contradicting the 
fact that u E M. □ 


5. FINITE ELEMENT FORMALISM 


We ntilize the variational principle nsed in Section 3 and a hnite element formalism to 
compnte the solntion pair (u, k) to the problem fll.Sp . for a prescribed energy flnx fll.6l) . 
This is essentially achieved by approximating the solntions to the constrained minimization 
problem (13.51) . 

Recall the admissible class A, dehned in (13.dh . Let R be a snbset of A, composed of N 
linearly independent fnnctions, Dehne the inner prodnct as 

rR 

{u,u) = 271 / ruudr, u,ueA, (5.1) 

Jo 

whose form is snggested by the constraint fnnctional fl3.3l) . Under the inner prodnct fl5.1l) . 
the set V can be orthonormalized via the Gram-Schmidt procednre. We let the fnnctions 
^ t)e orthonormal with respect to the inner prodnct fl5.ip . We approximate 
fnnctions u E Ahy nsing the hnite element formalism 

N 

“ = (5.2) 

i=i 


with oi,..., Oat G M. Using this formalism (15.2p . the constrained minimization problem (13.5p 
becomes 


min 


F{a) = X 



^ ^ oJj — Qo) Oj G 


(5.3) 


where a = (oi,..., oat), is called the variational vector. Note that X is a continnons, real- 
valned fnnction dehned over the snrface of the A-sphere of radins y/Qo centered at the 
origin in Hence, the constrained minimization problem fl5.3p is well-dehned and has a 
nontrivial solntion. 
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We use MATLAB’s Optimization Toolbox 2^ and the Chebfun package 1^ to solve 
fl5.3p . In particular, we obtain a minimum using the objective function as 

N 

2 


1 ^ fR\ ^2 


n" , , 2 a ^r'^'ipi'ipj^r 


i,j=i 


N 


1 + a X] 


k=l 


>dr 


Qo 


2 aTi 


In order to compute the wave propagation constant, we use a Lagrange multiplier A G M, 
such that {X{u),u) = X{{u),u). More explicitly, there exists a A G M such that 


-R 


n 


rUrUr H- UU 


2 ru^u 1 A \ ~j 

dr = An A / ruudr 


1 + au^ 


(6.4) 


JO L ' “r ) JO 

for every u G M. Comparing the weak formulation of the n-vortex equation (ll.bh . we get 
K = —27rA. Using u = u m. (15.dh and the prescribed energy flux Qo > 0, gives 

2rM"^ 

<5 ru; H- u" - - 

>0 

As in 


TT f 2 2 

K = -y:- / \rU^ + —U - -^ 

Qo An I r 1 + au^ 


dr. 


(5.5) 


36l |. we consider the case when the saturation constant, vortex winding number. 


and distance from the vortex core are: a = 0.1, n = 1, and R = 8, respectively. With this 
particular choice of parameter values and using Theorem 2.1 and Theorem 3.2{iii), we get 


0 < R < 9.9470. 


Equivalently, using fl5.5p . we get that the prescribed energy flux satisfies 


13.6 < Qo < cxD. 


(5.6) 


(5.7) 


The inequalities (15. 6 p and (15.7p are the necessary conditions for postive exponentially decay¬ 
ing solutions. This in turn, numerically demonstrates that Qo > |, as remarked in Section 
3, does imply that k > 0. 

Figure [H shows the soliton’s amplitude for several values of the energy flux Qq. Note 
that, as the prescribed energy flux Qq is increased the soliton’s amplitude also increases. 
The numerical error is estimated by substituting the formalism (15.2p into 


r-R 


n 


error = / I {ruQr - u + 2r 


-- — 2Kru dr. 

1 -I- / 


(5.8) 


For Qo = 40, 60, 80,100, we compute the propagation constant k and obtain n = 1.4901, 
2.5827, 3.2955, 3.8120 with error = 0.0001,0.0050,0.0120, 0.0116, respectively. Figure [H 
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FIG. 1. Soliton’s amplitude for a = 0.1, n = 1, i? = 8, and N = 40. 

also illustrates the behavior of the soliton’s amplitude for the borderline values of Qq = 
10,13.6, 20. The following values for the propagation constant k = —.0330, 0.0001, 0.0712 
with error = 0.0170, 0.0178, 0.0156, respectively, are obtained. As expected, the value of 
the wave propagation constant is negative when Qq < 13.6. 



FIG. 2. Propagation constant k as a function of the prescribed energy flux Qq with fixed parameter 
values a = 0.2, R = 8, and A = 15 for n = 1, 2,3. 

We also analyze the behavior of the solution pair {u, k) for a hxed value of the energy 
flux Qo by varying the vortex winding number n (see Figure [2]). Particularly, when a = 
0.1 and R = 8 , Theorem 2.1 states that the wave propagation constant k, must satisfy 
K < 10 — (rg + n^)/128, which imposes an upper limit on the vortex winding number of 
\n\ < v/1280 - rg ^ 35.6962 (i.e., | n\ < 35). However, for exponentially decaying solutions, 
i.e., K > 0, owing to Theorem 3.2{ii), the vortex winding number is bounded above by Qo/tt. 
Consequently and for example, using Qq = IOtt, we get |n| < 10 as a necessary condition for 
positive exponentially decaying solutions. 

Figure [3] shows the values for the wave propagation constant n = 0.7933, 0.0607, —0.4812, 
—0.8562, —1.3046 for n = 1, 2, 6, 8,10 with Qq = IOtt. We observe that the wave propagation 
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constant decreases as the vortex winding nnmber increases, which is expected and implied 
from the necessary condition of Theorem 2.1. In particnlarly, k —)■ —oo as n —)■ cxo. 

We remark that our numerical approach is in contrast with that of Skryabin and Firth 


36l |. We compute the wave propagation constant for a prescribed energy flux (see Figure 


On the other hand, Skryabin and Firth in [3^, compute the soliton’s amplitude for a 
prescribed propagation constant and then use fll.6l) to determine its corresponding energy 
flux. 



FIG. 3. Soliton’s amplitude for re = 1, 2,6, 8,10 with Qq = lOvr, a = 0.1, R = 8, and N = 20. 


6. CONCLUSION 


In this paper, we establish a series of existence results for ring-prohled localized opti¬ 
cal vortex solitons. We consider such solitons in the context of an electromagnetic wave 
propagating in a saturable nonlinear medium and model by a nonlinear Schrodinger equa¬ 
tion (II.ip . In particular, we focus on spatially localized ring-prohled optical vortex solitons 
governed by the re-vortex equation (II.bh . Below we summarize the results: 


1. From Theorem 2.1 and Theorem 3.2, a necessary condition for the existence of positive 
exponentially decaying solutions of the re-vortex equation (II.bh is 

re^ -|- Tq 


0 < K < a 


-1 


2i?2 


(6.1) 


Moreover, the vortex winding number must satisfy |re| < Qo/^r (see Theorem 3.2{ii)) 
and the prescribed energy hux Qq > 1/4 (see Theorem 3.1 («*)). Further, no small- 
energy-hux solutions exists for k > 0 when Qq < 1/4 (see Theorem 3.1 (m)). 


2. The existence of a positive solution is guaranteed by Theorem 3.1(i), however, the 
propagation constant k is undetermined. A lower bound for k is provided by Theorem 


2b 

















3.2(i), and an upper bound by Theorem 2.1. 


3. On a Nehari manifold, if the distance from the vortex core R is sufficiently large, 
then for any propagation constant satisfying fl6.1l) . a positive exponentially decaying 
solution exists (see Theorem 4.7 and Theorem 3.2{iii)). 

4. Using a finite element formalism, we compute the soliton’s amplitude and wave prop¬ 
agation constant for a prescribed energy flux. The numerical analysis shows that the 
wave propagation constant increases as the energy flux increases and decreases as the 
vortex winding number increases. Moreover, for given parameter values a, n, and R, 
we are able to numerically obtain a necessary condition for the existence of positive 
exponentially decaying solutions in terms of a prescribed energy flux (see fl5.7l) h 
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